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ABSTRACT 

TV-body simulations of Cold Dark Matter (CDM) have shown that, in this hierar- 
chical structure formation model, dark matter halo properties, such as the density 
profile, the phase-space density profile, the distribution of axial ratio, the distribution 
of spin parameter, and the distribution of internal specific angular momentum follow 
'universal' laws or distributions. Here we study the properties of the first generation 
of haloes in a Hot Dark Matter (HDM) dominated universe, as an example of halo 
formation through monolithic collapse. We find all these universalities to be present 
in this case also. Halo density profiles are very well fit by the Navarro et al (1997) 
profile over two orders of magnitude in mass. The concentration parameter depends 
on mass as c oc M°-^, reversing the dependence found in a hierarchical CDM universe. 
However, the concentration-formation time relation is similar in the two cases: earlier 
forming haloes tend to be more concentrated than their later forming counterparts. 
Halo formation histories are also characterized by two phases in the HDM case: an 
early phase of rapid accretion followed by slower growth. Furthermore, there is no 
significant difference between the HDM and CDM cases concerning the statistics of 
other halo properties: the phase-space density profile; the velocity anisotropy profile; 
the distribution of shape parameters; the distribution of spin parameter, and the dis- 
tribution of internal specific angular momentum are all similar in the two cases. Only 
substructure content differs dramatically. These results indicate that mergers do not 
play a pivotal role in establishing the universalities, thus contradicting models which 
explain them as consequences of mergers. 
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. 1 INTRODUCTION 
^ ■ 

5—1 ' mass distribution of the self-gravitating, quasi- 

d ' equilibrium dark haloes that form in an expanding uni- 
verse is an issue of fundamental importance. Early work 
on self-similar spherical collapse pre dicts virialized struc- 
tures with power-law den sity profiles ((Fillmore fc GoldreichI 
Il984l : [BertschingeilllQSSl ). As A'^-body techniques improved, 
it was realised that, in the hierarchical universes, the pro- 
file departs significantly from a single power law and is 
better fitted by a profile with curva ture in a log-log plot. 
lEfstathiou et al.| [l985: Dubinski fc Carlberd 11991. '). Later, 
it was established that the density profiles of haloes in CDM 
and other hierarchical clustering cosmologies, have a uni- 
versal for m which is well represented by the simple fitting 
formula of iNavarro et all l|l996l . Il997l . hereafter NFW) 

Ps 



where is a characteristic radius where the logarithmic 
density slope is —2, and Ps/4 is the density at r^. The 
above equation implies that p cx in the inner regions 
and p cx in the outskirts. A useful alternative param- 
eter for describing the shape of the profile is the concen- 
tration parameter c = r2oo/rs (r2oo is the virial radius de- 
fined as the radius within which the mean density is 200 
times the c ritical value) . Shortly after the NFW papers, 
IHuss et all (|l999bl ) showed that the NFW model also fit ha- 
los which form by monolithic collapse (i.e. without mergers) 
in a HDM cosmology. Nevertheless much work has focused 
on the role of mergers in establishing the NFW profile, since 
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during repeated mergers, concluding that the universal pro- 
file is generated by tidal stripping of small haloes as they 
merge into larger objects. iDekel et al.l l|2003l) extended this 
model and found the tidal compression from a halo core 
makes the satellite orbits decay from the radius where p cx r 
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to the halo center and causes a rapid steepening of the 
inner profile to p oc r~"[a > 1). Recent work has con- 
firmed the NFW hypothesis, extending the comparison to 
a much wider range of densities and radii than the orig- 
inal work, and uncovering small but significant deviation, 
between the mean densi ty profiles of simu l ated dark halos 
and t he NFW formula (IPower et al, 
20041: iMerritt et akl [ioo 



2003; Dicma nd et al 
Graham et al. 2006; Gao et al 



200a : lHavashi fc Whitell2007l ') . The differences, however. 



quite small compared to the scatter between different haloes 
of the same mass. 

In addition to the density profile, other halo proper- 
ties are found to follow universal profiles or universal dis- 
tributions in the hierarchical CDM model. The "phase- 
space density" p(r)/(T^(r) has a re markably a ccurate pure 
power-law distribu tion with radius (i Tavlor fcN avarro 200ll; 
iBarnes et al.|[200^ . The distributions of shape parameters 
(e.g. the axis ratios) have a weak dependence on mass and 
redshift llBuUocldl2002|. iKasun fc Evrardl2005l : lAlkood et aD 



l2006l : iBett et al.ll2007l ). The spin parameter distribution is 
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mass distribution of specific angular momentum j is well fit- 
ted by a universal function. All these distributions appear 
to depend little, if at all, on the global cosmological parame- 
ters and on the shape of the initial matter power spectrum. 
However, the origin of these universalities is still not well 
understood. 

In this paper, we will test whether mergers are the dom- 
inant physical mechanism to produce these universal profiles 
or distributions. To this end, we study a range of properties 
of the haloes which grow by two different paradigms, by 
hierarchical aggregation and monolithic collapse. The con- 
cordance cosmology (ACDM) is a standard hierarchical uni- 
verse, and the HDM dominated universe provides an ex- 
ample where formation of the first haloes is monolithic. In 
a HDM universe, small objects cannot form because free- 
streaming effects smooth out small-scale structure in the 
initial condition. The first halos then form by smooth col- 
lapse. 

We begin in section [2] with a summary of the A^-body 
simulations and halo catalogues used in this paper. In sec- 
tion [3l we give a general description of how halos form and 
evolve in the HDM universe. In section U) a range of halo 
properties are studied in depth in our two model universes. 
Then, we will summarise and discuss our results in section O 



2 THE SIMULATIONS AND HALO 
CATALOGUES 

In this work, we simulate the HDM density field using 
512"^ particles within a 100h~^Mpc cube. For simplicity, we 
choose an Einstein-de Sitter universe dominated by a sin- 
gle massive neutrino. Then the cosmological matter density 
parameter is Q,m = 1 and the simulation particle mass is 
2.07 X lif Mq/H. The initial power spectrum used to per- 
turb the i nitial particle se t is ba sed on the theoretical pre- 
diction of iBond fc Szalavl ((198^^). The power spectrum in 
a HDM universe possesses a coherent free streaming scale: 




Figure 1. The evolution of a massive halo in a HDM universe 
with mass 1.3 X lO^^MQ/h at redshift zero. All particles within 
^"200 at 2 = are traced back to seven higher redshifts. In addition 
the particles in the very inner region r < 0.01r200 a^re traced 
back and are shown as the blue points. The size of each box is 
8Mpc/h in physical (not comoving) units. The redshift and the 
corresponding eg are showed in each panel. 



5-nvL ~ 22.2Mpc in our case. The corresponding damping 
mass IS 7.28 X 10^"* Mq/Zi. This power spectrum is normalised 
to CT8 = 2 which corresponds to the collapse of the first non- 
linear structu res at z ^ 6. This sim ulation was presented 
as glass512 in IWang fc Whiti (|2007l . hereafter paper I) and 
more technical details can be found there. We choose to 
compare it with a simulation of the WMAPl model from 
IWang et all (|2008l '). This simulates a concordance ACDM 
universe with 540'^ particles in a 125h~^Mpc cube. Hereafter 
we refer these two simulations as HDMRUN and CDMRUN. 
The softening lengths for these two simulations are 10 and 
5kpc/h in the HDM and CDM cases respectively. 

As pointed out in Paper I, the population of haloes 
in HDMRUN is contaminated because discreteness effects 
break up the filaments and produce many artificial small 
haloes. The resolution limit above which haloes are immune 
to this effect is Alum ~ 8.8 x W^^Mq/H, and corresponds 
to about 5000 particles. In the current study we therefore 
focus primarily on well resolved haloes with particle number 
A*'2oo > 10000 within virial radius r2oo in two simulations. 
Several of our most massive HDM haloes experience one 
or two major mergers at low redshift as in the CDM case. 
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Figure 2. Mass assembly histories (MAH) for the 10 most mas- 
sive haloes at z = in the HDMRUN sample. The red curve 
indicates the halo presented in Fig. [T] 



We exclude these second generation haloes from our anal- 
ysis below. The number of haloes which satisfy the above 
constraints is 58 and 304 in HDMRUN and CDMRUN re- 
spectively. When studying the Af200-c relation, we reduce 
the particle number limit to 5000 in order to cover a wider 
mass range. We then have 84 (HDMRUN) and 1752 (CDM- 
RUN) haloes in the two simulations. 

Our halos are identified b y a standard b = 0.2 friends- 
of- friends (FOF) group-finder (Davis ct ahl ligSSl ). Then the 
SUBFIND algorithm (Springcl ct al. 2001) is used to resolve 
these objects into substructures and main subhaloes. The 
latter then define our halo sample through our A^200 lim- 
its. At high redshift, the FOF method tends to link more 
particles together especially in our HDMRUN where large 
filaments and sheets are easily joined. SUBFIND defines the 
centre as the minimum of the gravitational potential, and 
this is used to estimate A^2oo ■ We also use merger trees built 
up from the s ubhalo catalogue pro duced by the SUBFIND 
algorithm fsee ISpringel et aLlbOOSl ). 



3 MONOLITHIC GROWTH 

In Fig: [1] we study the evolution of a typical halo in our 
HDMRUN simulation. The mass of this halo is about twice 
the free-streaming mass scale. In the lower right corner 
[z = 0), all particles within r2oo are chosen. These are 
then traced back to the initial condition in the other panels. 
The particles within the central region r < 0.01r2oo are also 
traced back to high redshift and are highlighted with blue 
points. It is interesting that these images are not similar to 
those found in the con cordance CDM u niverse using a sim- 
ilar representation bv Gao et al.l (|2004h : the whole density 
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Figure 3. Formation histories of two "first" objects in a HDM 
universe. The right column shows a halo with N200 = 2686 at 
z = 2.9 (roughly M200 = 3.5 X lO^M©); while the left column 
shows another halo with A''200 = 1215 (M200 = 1-6 x 10^"' Mq) at 
z = 4.3. All particles in both halos at the last redshift are traced 
back to higher redshifts and presented as blue points. Particles in 
the very inner region r < 0.1r200 are traced back and presented 
as red points. Other background particles are shown with black 
points. The coordinate is comoving scale, and the size of each box 
is 6Mpc/h. The redshift is labeled in each panel. The simulation 
glassl28 (see Paper I) is used here. Note the artificial discreteness 
features visible in the filaments in the z = 3.6 panel. 



field is smooth at all times and no obvious substructures are 
seen. The collapse along the filaments into knots occurs si- 
multaneously with the accretion of diffuse particles. At red- 
shift 1.14 and 0.84, some parts of the filaments are missing 
since they fall outside r2oo at 2 = 0. The particles which end 
up within 0.01r2oo stayed close together at all times, even 
in the initial condition. They fall to the centre in a smooth 
sp herica l collapse. In a CDM universe, as presented in Fig 2. 
of Gao et al. (2004), the matter which ends up in the central 
region typically comes from a number of different objects at 
early times. 

In Fig[5]we present the mass assembly history (MAH) of 
the 10 most massive haloes in the HDMRUN sample. This is 
just the mass growth of the most massive progenitor in each 
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case. We find that the MAHs are similar to those of haloes 
in a CDM universe: at early stages, the mass grows rapidly; 
but this slows dramatically towards 2 = 0. For example, the 
MAH of the halo shown in Fig. [T] is presented as a red curve 
in Fig (2] This halo stays in the rapid growth phase until z — 
1.2, and its assembly is dominated by smooth spherical infall 
as Fig. [T] indicates. This appea rs quite different from the 
hierarchical growt h described bv lWechsler et ahl (|2002l ) and 
IZhao et al.l l|2003l ) their rapid growth phase is dominated 
by mergers. They also speculated that the universal inner 
density profile results from violent relaxation during this fast 
merger phase. We will see below that the density profile is 
also universal in our HDM simulation. Thus the NFW profile 
apparently does no t require mergers, as noted originally by 
iHuss et al.l (Il999al l. 

Halo formation histories obviously differ in CDM and 
HDM universes. Another interesting question is how and 
where the first objects form in the HDM case. Here we must 
keep in mind that the free-streaming scale in HDM is just 
a characteristic scale and its relationship to the actual mass 
of the first objects is not obvious. In Fig. |31 we present for- 
mation histories for two "first" objects. They were identified 
at z = 4.3 (left column) and z = 2.9 (right column) when 
they had gathered 1215 and 2686 particles within r2oo re- 
spectively. In order to see the particles better, we here use 
a lower mass resolution simulation glassl28 which was also 
presented in Paper I. This simulation has the same initial 
density field as HDMRUN, but 16 times lower mass resolu- 
tion. As in Fig [T] we trace back all particles within r2oo 
(blue) and within 0.1r2oo (red) to higher redshift. Other 
background particles are shown as black points. It is obvious 
that both haloes form by smooth spherical collapse. In other 
words, the particles in the filaments and "voids" fall into the 
halo together. Even the particles within the central region 
(< 0.1r2oo) seem to follow a roughly spherical collapse. This 
agrees with what we see in Fig[T] In the right column, the 
halo forms in a node which connects three filaments, and the 
halo forms almost at the same time as the three neighbour 
filaments. For the case in the left column, the halo forms 
at the end of one filament {z = 6.0) and then grows at the 
same time as another filament. We have also checked some 
other examples and find that almost all of them are born 
in nodes or at the ends of filaments and grow by spherical 
accretion. 

Merger events in such a monolithic universe are ex- 
pected to be rare. We find this to be true in our HDM 
simulations. For example, in HDMRUN, if we define merger 
events as the merging of 'real' haloes whose masses are larger 
than the mass limit, only a few haloes experience a major 
merger event (with mass ratio of the two progenitors greater 
than 1:3). We exclude these haloes from our study in this 
paper. Only about 30 percent of the remaining haloes ex- 
perience minor merger events, and it is thus reasonable to 
regard them as a sample of 'first' haloes. 



4 HALO PROPERTIES IN THE TWO 
UNIVERSES 

In this section, we will study a range of "first" halo prop- 
erties in the HDM universe and compare them with their 
counterparts in the ACDM model. 
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Figure 5. The distributions of halo mass M200 and concentration 
parameter c. The top panel is from HDMRUN and the bottom 
one from CDMRUN. 
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Figure 6. The distribution of formation time for haloes with 
M200 ^ 10^3 M0/h in the HDMRUN sample. 



4.1 Density Profile and Formation History 

Previous studies jMoore et al.l 1 19991 : IColfn et al] |200G| : 
lEke et al.ir200ll : iBusha et al.ll2007l ) have focused on the for- 
mation of the first haloes in a WDM cosmology and have 
claimed that the density profiles of these haloes do not dif- 
fer substantially from the NFW form found in a CDM uni- 
verse. We support these findings also in our HDM case, con- 
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Figure 4. Density profiles for 16 haloes with particle number within r2oo (^20o) greater than 10,000. These 16 haloes cover two orders 
of magnitude in mass. The density profiles are normalized by r'^/pcrit ^'Hd the radius are normalized by r2oo- Iri each panel, the red 
dotted curve is the NFW fit to the numerical measurement (black solid curve). Two vertical dotted lines show the softening length and 
r200- The corresponding M200 (in unit of h~^MQ ) and the concentration parameter c are listed in each panel. 



firming the original result of lHuss et al.l (|l999lj ). In Fig. |4j 

we present density profiles for 16 haloes with mass from 
3 X IO^^Mq/Zi to 4.5 X IO^^Mq/Zi. All these haloes include 
more than 10, 000 particles. Here the concentration param- 
eter is measured by fitting an NFW profile to the numerical 
results using logarithmically spaced radial bins in the range 
2e < r < r2oo. The softening length e and r2oo are shown 
in the plot by vertical dotted lines. We find these profiles to 
follow the NFW model very well. It is also obvious that the 
ps and c have a strong mass dependence. Both parameters 
increase with increasing halo mass. These mass dependences 
disagree with those in a CDM universe where more massive 
halos have a lower ps and c. In Fig. (5] we display the de- 
pendence of concentration parameter on halo mass for our 
two halo samples. In the lower panel, we present results for 
CDMRUN. Th e mass dep e ndenc e, c o c M^n'^^, is c l ose to 
that found by iNeto et al'l l|2007f) and iMaccio et ahl (|2007l ') 
c oc M^q'^^ . In the upper panel, we present this relation 
for HDMRUN. It is interesting that the mass dependence 
inverts and follows c oc M200: the more massive a halo, the 
larger its concentration parameter. 

Many previous studies have found that the structural 
properties and the mass accretion his tories of haloes are 
closely related in a CDM universe (e.g. N avarro et al.|[T996l . 
[1997, ; Wechsler et a l. 2002; Zhao et al. 2003|). Concentration 
increases with the formation time and the characteristic den- 
sity ps can be related with the mean cosmic density at the 
time of formation. In Fig. [5] , we check this relation for our 
HDMRUN haloes. The formation time is defined here as the 
earliest time when half of the halo mass was in its main 



progenitor. It is interesting that massive haloes form a bit 
earlier than their low mass counterparts in this cosmology. 
Combining with the M200 ~ c relation showed in Fig. [S] we 
find that the earlier forming objects do indeed have a larger 
concentration parameter, in agreement with the result for a 
CDM universe. This indicates that here also the inner part 
is assembled during the fast early growth phase. 

4.2 Kinematics 

If haloes follow a universal density profile as discussed 
above, and are in an equilibrium state, then the solu- 
tion of Jeans equation for spherical, isotropic systems indi- 
cates that their kinemat ic structure may also be universal. 
I Taylor fc Navarrol (|200ll 'l have found that, for dark haloes 
in a CDM universe, there is a power-law relationship be- 
tween "phase-space density" and radius: p/o""^ ~ r~°' with 
a = 1.875, where phase-space density is defined as the ratio 
of local matter density p to the cube of the local veloc- 
ity dispersion a. This phase-space density is inversely re- 
lated to the local entropy density. In the semi-analytic ex- 
tended secondary infall model, this nearly scale-free nature 
of p/cr^ is a robust feature of virialized haloes in equilibrium 
I Austin et al.ll2005h . Further investigation of halo formation 
processes indicates that this scale-free feature cannot be the 
result of hierarchical merging; rather it must be an outcom e 
of violent relaxation (I Austin et al.|[2005l : lBarnes et al.ll2006l ). 
We show results for massive haloes in our HDM and CDM 
samples in Fig. [7] In order to reduce the noise, we stack the 
20 most massive haloes in each case. Before the stacking. 
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Figure 7. Stacked phase density profiles for tlie 20 most massive 
haloes in our two halo samples. The profile of each halo is normal- 
ized by a factor ffoo/Pc^t before stacking. The dotted straight 
lines are our power law fits with indices a = 1.94 and a = 1.83 
for HDMRUN(red triangles) and CDMRUN (blue squares) re- 
spectively. In order to get a better dynamic range, all symbols 
and lines multiplied by {r/r2oo)^- The error-bars indicate the la 
scatter. 



the profile for each halo is normalized by a factor Vioo/ pcrit- 
The error-bars in the plot indicate the Icr scatter. We find 
the two profiles to follow a power laws moderately well. The 
fitted indices differ slightly: a = 1.94 and 1.83 for HDMRUN 
(red) and CDMRUN (blue) respectively. 

In addition, the velocity anisotropy in haloes depends 
on radius: The dispersion tensor is isotropic near the center 
and moderately radially anisotropic near the virial radius. 
In Fig. |8l we present averaged velocity anisotropy profiles 
/3(r) = 1 — <Jt/(^r for the 20 most massive haloes in our 
CDM and HDM samples. Here at and ar are the tangential 
and radial velocity dispersions respectively. We again find 
very similar results for the two cases. Outside 0.2r2oo, the 
mean radial anisotropy in HDMRUN (red) is slightly larger 
than in CDMRUN (blue), but the effect is very small. 

The similarities in phase-space density profile and ve- 
locity anisotropy profile in the two cases indicate that these 
properties are also universal and depend little on whether a 
halo is assembled by mergers or by monolithic collapse. 



Figure 8. The stacked velocity anisotropy profile /3(r) for the 20 
most massive haloes in the HDMRUN (red triangles) and CDM- 
RUN (blue squares) samples. The error-bars indicate the Icr scat- 
ter. 
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4.3 Halo Shape 

In Fig. |9] we show axis ratios for haloes in HDMRUN and 
CDMRUN. We define the axes using the inertia tensor of 
the mass distribution within the virial radius r2oo- 



lal3 



1 



(2) 



Figure 9. Scatter plots of the axis ratios of haloes in HDM- 
RUN and CDMRUN (top two panels) and histograms of their 
probability distributions (lower two panels: red dotted and black 
solid curves are histograms for CDMRUN and HDMRUN respec- 
tively). Here a, b and c are the major, intermediate and minor 
axis lengths. Kolmogorov-Smirnov tests indicate that the distri- 
butions in two universes do not differ significantly. 
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Figure 10. The probability distribution of the spin parameter 
A' in HDMRUN (blue solid) and CDMRUN (red dotted). The 
black solid curve is the fit to a CDM universe from lBuUock et al.l 
1I2OOII) . 



Figure 11. Mass distributions of specific angular momentum 
within 58 massive haloes in HDMRUN. Each coloured symbol 
presents one halo. Each mass distribution is normalized by its 
own virial mass (Mq) and shape parameter fi, all symbols then 
stay around the black curve jtt— • 



where the ri are the positions of the Np particles within 
r2oo- a and /3 are tensor indexes with values of 1,2 or 3 and 
indicate the three components of each particle's position. Af- 
ter diagonalizing this matrix, characteristic axis lengths are 
found as the square root of the eigenvalues. In the top two 
panels of Fig. |9l the intermediate-major (b/a) and minor- 
intermediate (c/b) axis ratios of halo samples from the two 
simulations are presented. The distributions of shape pa- 
rameters are clearly very similar in the two cases. The mean 
values of the axis ratios b/a and c/b are both close to 0.82. In 
the lower two panels, we compare the probability distribu- 
tions of b/a and c/b in the two simulations. The Kolmogorov- 
Smirnov test shows no evidence for any difference between 
the distributions in the HDM and CDM cases. 

The above results indicate that different growth paths 
(monolithic or hierarchical) affect the global shape of halos 
very little. 



4.4 Spin and Angular Momentum Distribution 
within Haloes 

The conventional measure of halo angular momentum, the 
dimensionless spin parameter A, is defined in terms of mass, 
energy and angular momentum and is related to the ratio 
between a halo's mean angular velocity {u>) and the angular 
velocity w hich would be required to support it by rotation 
alone (cjo) (|Peebleslll969l ): 



A ; 



J\E\'/^ 



4— 



(3) 



The total angular momentum J and energy E are needed to 
calculate this parameter. However, in this study, we follow 
iBuUock et all l|200ll ) and define a more easily measured spin 
parameter A' as: 




10.00 



Figure 12. The probability distribution of the shape parameter 
^l in HDMRUN (blue solid) and CDMRUN (red dotted). The 
black solid curve is the fit to a CDM universe from lBuUock et al.l 
l|200ll 'l. 



A' = 



J 



\/2M2ooVr2oo 



(4) 



where J is the total angular mo mentum of all particles 
within r2oo and V = y^?M20o/»^oo is the circular velocity at 
radius r2oo. This definition gives similar value to Equ.|3] Be- 
cause A' is defined for isolated systems but applied to haloes 
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in their cosmic context, the practical definition of a halo is 
more critical than the choice of definition of A l|Bett et alj 

The nett spin of a dark halo is acquired fro m torques ex- 
erted by neighbouring structures at early time s (|Hoylell 19491 : 
iPeeblesll 19691 : lDoroshkevi"cblll970l : Iwhitell 1984 ) . and does not 
evolve much after the turnaround point of a halo's MAH. Af- 
ter this time, the moment of inertia of the collapsing material 
decreases and the universal expans io n reduces the strength 
of tidal forces l|Porciani et al.ll2002l ). IVitvitska eTe^ (|2002l ') 
noted that the spin parameter fluctuates strongly with time, 
depending on the details of assembly: the spin increases 
abruptly during a m ajor merger and decreases gr adually be- 
tween such mergers. iD'Onghia fc Navarro! l|2007h found this 
effect to reflect the unrelaxed nature of the system; equi- 
librium haloes show no significant correlation between spin 
and merging history. In our IfDM simulation, most haloes 
grow in a monolithic way, so the impact from major mergers 
should be negligible. 

In Fig. 1101 we compare the distribution of spin param- 
eters in our HDMRUN halo sample with that in the CDM- 
RUN sample. A Kolmogorov-Smirnov test shoes no signifi- 
cant diffierenc£_between_two different cases. 

iBullock et al] l|200lD noticed that the distribution of 
specific angular momentum j within the CDM haloes has 
a universal profile, specifically that the distribution of mass 
over specific angular momentum j is well fit in most haloes 
by the universal form: 



100.0 



M(< j) = Mo 



JO + j 



(5) 



where Mo is the virial mass, which we here replace by M200. 
jo = {11— l)jmax and jmax IS the maximum specific angular 
momentum. This behaviour is not truly universal since fi 
is an adjustable shape parameter which varies from halo to 
halo. We now check this distribution for haloes in our HDM 
and CDM simulations. 

In order to resolve the j profile adequately, we only 
consider haloes with particle number larger than 3 x 10'*. 
Different regions are defined as cells in the usual spherical 
coordinates {r,9,(f>). We make sure that each cell contains 
an approximately equal number of particles. We first di- 
vide the whole halo into 10 radial shells which host almost 
equal particle numbers. Then each shell is further divided 
into six zones: each cell spans the full 2-k range in (j> and 
spans the equal solid angle between cos(S) = —1 and 1. 
The two zones with the same r and |cos(6')| that are above 
and below the equatorial plane are assigned to one cell. In 
Fig. 111! we present the distribution of normalized mass frac- 
tion M(< j)/(Mo/i) for 58 haloes from HDMRUN. Different 
coloured symbols present different haloes. We find that, after 
the parameter fi has been adjusted, all haloes are close to the 
curve x/{l + x), thus the specific angular momentum profile 
can be described by one shape parameter /x very well. We 
compare the probability distribution of this parameter for 
the HDMRUN and CDMRUN samples in Fig. [12] The dis- 
tributions a re similar and c lose t o a log-normal distribution 
as noted bv lBuflock et al.l (|200lD ( A Kolmo gorov-Smirnov 
test show s no indication of a d ifference.) A similar result was 
found by IChen fc Jin"3 l|2002f ) for WDM haloes: the distri- 
butio n is statistically indist inguishable from the CDM case . 
As in Ivan den Bosch et al.l (12002 ) and lChen fc Jind l|2002h . 



10.0 



A 




10"^ 10""^ 10"^ 10" 

Msub/M200 



Figure 13. The stacked cumulative subhalo mass function of the 
30 most massive haloes for our two samples: red diamonds and 
blue stars are for HDM and CDM haloes respectively. The results 
for several individual HDM haloes are also presented as black 
curves. The HDM curves are divided into two parts at 0.1M;j^; 
real (spurious) subhaloes have masses greater (less) than this 
mass limit and are shown with solid (dotted) curves. The ver- 
tical dashed line indicates the mean mass limit for all 30 haloes 
in HDMRUN < O.IMh™ >. 



we also find some cells with negative angular momentum. 
These are excluded from our analysis. 

The similarity between the distributions of spin pa- 
rameter and internal specific angular momentum in the 
HDM and CDM cases shows that formation through merg- 
ers is not necessary to generate "universal" angular momen- 
tum distributions. Th is confiicts with the explanation which 
IVitvitska et al. I liooa ) give for the origin of the angular mo- 
mentum distribution within CDM haloes. 



4.5 Substructure 

In Paper I, we gave a detailed discussion of the formation 
of substructures in a HDM Universe. We found almost all 
substructures in HDM haloes are either a result of spurious 
filament fragmentation (at low masses) or major mergers at 
late times (at large masses). As we see in Fig[TJ ai z — 0.84 
and z = 0.58, there are many artificial regularly spaced ob- 
jects falling into the clusters from filaments and these be- 
come substructures of the HDM haloes. This effect prevents 
us from studying the abundance of the real subhaloes at low 
mass since it is hard to identify which subhaloes (if any) 
are real. In Fig |13l we present cumulative mass functions of 
subhaloes for stacks of the 30 most massive haloes in each 
sample. For the HDM case, the subhalo mass functions of a 
few individual haloes are presented with black curves. For 
haloes in the HDM universe, we define 'real' subhaloes to 
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have mass greater than O-lAfum = 8.8 x 10^^ Mq /h, and 
smaller subhaloes are considered 'spurious'. The mass func- 
tion of the 'real' and 'spurious' subhaloes are shown as solid 
and dotted lines respectively. The vertical dashed line indi- 
cates the mean mass limit < Q.lMum. > for all 30 haloes 
in the HDMRUN. We have assumed here that the average 
mass loss fraction is 0.1 when a subhalo falls into its parent 
halo and suffers tidal stripping. It is interesting that both 
distributions follow a similar power law. That of the HDM- 
RUN sample is approximately one order of magnitude lower 
than the counterpart in CDMRUN and on average very few 
'real' subhaloes are found in each halo in the HDM case. 
This reflects the lack of hierarchical build up in the HDM 
case. 



5 CONCLUSIONS AND DISCUSSION 

In this paper , we study the haloes which form by monolithic 
collapse in HDM cosmological simulations, and we compare 
a range of properties of these haloes with these of their coun- 
terparts in a ACDM hierarchical universe. From this com- 
parison we explore which physical mechanisms are respon- 
sible for the universal profiles or distributions of halo prop- 
erties. Our simulated HDM universe has an inherent char- 
acteristic scale below which the formation of small haloes 
is suppressed. In this universe, the first generation of halo 
forms in a "top-down" way - more massive haloes form ear- 
lier by smooth accretion. Mergers occur in significant num- 
bers only when building up later halo generations. The HDM 
cosmology thus provides a good example of monolithic halo 
growth from Gaussian initial condition. Furthermore, be- 
cause the HDM characteristic scale is much more sharply 
imprinted than in the WDM case, we have a large and well 
defined monolithic halo sample. As a result, we can explore 
monolithic growth in a better controlled way than previous 
work based on WDM universes. We summarise our conclu- 
sions on monolithic halo formation as follows: 

(i) First generation haloes grow by smooth accretion. 
Such haloes form from roughly spherical collapse of matter 
around the nodes or ends of filaments. The formation his- 
tories of these monolithic haloes are also characterised by 
two phases: fast initial and slow later accretion, very sim- 
ilar to those found in a CDM universe. In the fast growth 
phase, smooth infall is the main mechanism rather than the 
mergers seen in the CDM case. 

(ii) The density profiles of our monolithic haloes are 
well fit by the NFW profile. The concentration parameter 
c and the characteristic density ps have strong dependences 
on halo mass: the more massive a halo, the larger c and 
Ps- This is the inverse of the CDM dependence, but the 
concentration-formation time relation is quite similar in the 
two cases: earlier formed haloes tend to be more concen- 
trated than their later formed counterparts. 

(iii) Phase-space density profiles and velocity dispersion 
anisotropy profiles are very similar in our two cases. This 
indicates that the kinematic structure of haloes is generated 
by wide variety of dynamical collapse processes, not just by 
mergers. 

(iv) The distribution of shape parameter is very similar 
in HDM and CDM universes, as are the distributions of spin 
parameter A' and of internal specific angular momentum. 



(v) In our HDM simulations, most subhaloes result from 
the infall of spurious small haloes which form in the filaments 
as a result of numerical discreteness. The subhalo abundance 
is much lower than in CDM haloes, and indeed there are very 
likely no subhaloes in first generation HDM haloes. 

The above results show that, except for substructure, 
the properties of haloes formed by monolithic collapse are 
very similar to those of haloes formed by hierarchical cluster- 
ing in a CDM universe. In particular, they have very similar 
profiles of density, phase-space density and specific angu- 
lar momentum, and very similar distributions of axial ratios 
and spin parameter. These results indicate that mergers are 
not responsible (or are not required) for the universal prop- 
erties of CDM haloes. This agrees with results from earlier 
work based on constr ained simulations (|Huss et al ] |l999al : 
iMacMillan et al1l2006l ). Attempts to explain these universal 
properties with merger-driven models seem unlikely to be 
on the right track. 

The results in this study could also find some appli- 
cations in realistic dark matter models, especially for the 
WDM cosmology. In our HDM universe, first generation 
haloes form by smooth near-spherical accretion with ini- 
tial mass well below the characteristic free-streaming mass. 
These haloes form simultaneously with the filaments and 
sheets. This picture is different from that in the familiar 
CDM universe. Our results suggest that the structure of the 
first haloes, in a WDM universe (and thus, presumably, also 
of their later descendants) will differ little from the structure 
of CDM haloes, despite the differences in assembly history. 
This makes it seem unlikely that changing to a WDM model 
will help resolve the apparent problems of the CDM model 
on galactic scales. 
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